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Using molecular dynamics simulations, we investigate systematically the water permeation prop¬ 
erties across the single-walled carbon nanotube (SWCNT) in the presence of the terahertz electric 
field (TEF). With the TEF normal to the nanotube, the fracture of the hydrogen bonds results 
in the giant peak of net fluxes across the SWCNT with a three-fold enhancement centered around 
14THz. The phenomenon is attributed to the resonant mechanisms, characterized by librational, 
rotational, and rotation-induced responses of in-tube polar water molecules to the TEF. For the 
TEF along the symmetry axis of the nanotube, the vortical modes for resonances and consequently 
the enhancement of net fluxes are greatly suppressed by the alignment of polar water along the 
symmetry axis, which characterizes the quasi one-dimensional feature of the SWCNT nicely. The 


resonances of water molecules in the TEF can 
designs used for various purposes. 

I. INTRODUCTION 

The study of the water molecular transport in 
nanochannels is important not only for the basic research 
and potential applications but also for the understand¬ 
ing of the water permeation in biological channels. In¬ 
deed, one important reason why research on nanochan¬ 
nel transport properties has continued to be of interest 
is that many primary characteristics of water confined 
in simple nanochannels are similar to that of complex 
biological channelsIn the past, the water trans¬ 
port properties have been studied extensively in single- 
walled carbon nanotubes (SWCNTs)®^^^. Novel studies 
on water transports in nanochannels have continued to 
appear in recent years with various designs concerning 
temperature gradients^^, charge modification^^, coulomb 
dragging^®, collective mode oscillations^®, and phonon- 
induced frictions^^. Besides the biological implications 
in these studies, the stream of these novel progresses 
may potentially become integrated into the application 
pool for the design of innovative nanofluidic devices, such 
as flow sensors^®, desalination of seawater^®, molecular 
sieves^®, and so on. 

Recently, broad interest has been aroused to study 
the electric field (EF) effects on the water permeation 
It is known that the applied EF can stem from 
electromagnetic radiations. A well-known example is the 
laser-generated EF that has fascinating applications, e.g., 
see®^’®®. Nowadays, high-power terahertz (THz) radia¬ 
tion can be generated progressively, say, by optical rec¬ 
tification methods®'’’, linear accelerators®®’®®, and laser- 
driven ion or electron accelerators®^’®®. The consequent 
THz EF (TEF) of the order of GV/m, generated by 
focusing high-power THz rays, can open up a new era 
to study structural transitions in polar molecules given 
the TEF-driven molecular rotations over random ther¬ 
mal motions. Such studies are of fundamental impor- 
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have potential applications in the high-flux device 


tance in revealing the mechanism for molecular motors 
that produce directional motion®®^'’^. While for water 
permeation in nanochannels, the study of TEF effects 
is scarcely found in the literature, and no attention has 
been paid to the role of the TEF-driven molecular ro¬ 
tations. It is thus of prime significance for us to first 
reveal whether and how the TEF-driven molecular ro¬ 
tations and intrinsic mechanisms affect the water flow. 
A clear understanding of the TEF effects is also very 
instructive to develop new devices for efficient water fil¬ 
tration and energy transfer. Using molecular dynam¬ 
ics (MD) simulations, we have observed that water per¬ 
meation properties can greatly be affected by the TEF 
through rotational and rotational-induced translational 
resonances in a broad THz frequency range. 


II. COMPUTATIONAL METHODS 

The MD simulation system built with the molecular 
visualization program^® is illustrated in Fig. 1. An un¬ 
capped (6,6) armchair SWCNT with a length of 1.34 nm 
and a diameter of 0.81 nm is embedded in two parallel 
graphite sheets along the z direction. The distance be¬ 
tween the bottom end of the SWCNT and the graphite 

O 

sheet is 2 A. MD simulations are performed using the 
NAMD2^^ in a NVT ensemble with the initial box sizes 
of Lx = 3.5 nm, Ly = 3.5 nm, Lz = 6.3 nm, con¬ 
stant temperature (300 K) achieved by the Langevin 
dynamics with a damping coefficient of 1 and pe¬ 

riodic boundary conditions in all directions. The elec¬ 
trostatic interactions are handled by the particle mesh 
Ewald method^® with a real space cutoff of 1.2 nm, and 
the cutoff for the van der Waals interaction is also 1.2 nm. 
The CHARMM27 force field^®’^^, TIP3P water model^® 
and SHAKE algorithm for bond lengths^® are adopted 
in all the simulations. The time step is 1 fs. To obtain 
a directed flow'’^, an acceleration of 0.01 nm ps~^ per 
water molecule is imposed along +z direction to produce 
a pressure difference between two ends of the SWCNT 
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about 20 MPa^'^. The carbon atoms at the inlet and 
outlet are fixed and other atoms of the SWCNT are flex¬ 
ible to prevent possible errors from the rigid wall^°’^°. 



FIG. 1: (Color online) Snapshot of the simulation 
system. An uncapped (6,6) armchair SWCNT 
combined with two graphite sheets is solvated in a 
water box of 3.5x3.5x6.3nm^ with 1852 molecules. 
There is a vacuum between the two sheets. The axis of 
SWCNT is parallel to the z-axis. To obtain the net flux, 
an acceleration of each water molecule is imposed, see 
text. 


In our simulations, only the space between two reser¬ 
voirs is exposed to the spatially uniform TEF, and a jus¬ 
tification for this setup is given in the Supporting Mate¬ 
rial®^. With the EF direction being parallel to the x-axis, 
the TEF, acted on the atomic partial charges, can be ex¬ 
pressed as 

E = £'o(cos(wt), 0,0). (1) 

In this work, we take the amplitude of the EF Eq = 
1 V/nrn^'^’^^ along the a;-axis, unless otherwise denoted, 
and consider the frequencies ranging from 0.1 to 40 THz. 
Hereafter, these far infrared EFs are still dubbed the 
TEF, since the interesting results lie below the region 
over a dozen THz. It is known that the electric and 
magnetic fields of the electromagnetic wave obey the re¬ 
lation |E/B| = c where c is the speed of light. Because 
the speed of water molecules inside nanotubes is far less 
than the speed of light, the Lorentz force is negligibly 
small in comparison to the electric force. Thus, we ne¬ 
glect the effect of the magnetic field in our estimation. 

We should point out that the fast water transport de¬ 
sired in this work is different from the in-tube water 
pumping by the AC electric field^^’^^. For water pump¬ 
ing through SWCNTs, the energy acquisition in the AC 
field is necessary, and simultaneously the spatial asym¬ 
metry should be temporally retained in a delicate ther¬ 
modynamic relaxation process. Our simulations with the 
symmetric setup in the uniform TEF are carefully testi¬ 
fied to have null pumping under zero pressure gradient. 


III. RESULTS AND DISCUSSIONS 

Simulations are performed for the system within the 
EF at different frequencies. For each frequency, the sim¬ 
ulation time is 45 ns, and the last 40 ns is collected 
for analysis. For clarity, here the net flux is defined 
to be the difference between the water molecular num¬ 
ber per nanosecond across the SWCNT from one end to 
the other^do. The average occupancy hydrogen 

bond numbers of water molecules inside the SWCNT are 
denoted by the symbol N and Nh, respectively. 

The net water fluxes, as a function of the EF frequency, 
are shown in Fig. 2 where the zero-field {Eq = 0) net flux 
(about 21 ns~^) is also displayed. We see that with the 
EF along the cc-axis the net fluxes can be characterized 
by a giant peak centered at 14 THz with a three-fold en¬ 
hancement of the flux. In sharp comparison to the zero- 
field result, the dramatic enhancement of the water flux 
is observed in the frequency range around 10 THz to 20 
THz. For much lower or higher frequencies, the net flux 
falls off, being close to the one of the zero field. For the 
net fluxes with the EF along the z-axis, also displayed in 
Fig. 2, the enhancement is greatly suppressed. This dis¬ 
tinct phenomenon, associated with the z-direction polar¬ 
ization of in-tube water molecules, characterizes interest¬ 
ingly the property of the quasi one-dimensional SWCNT. 



Frequency(THz) 


FIG. 2: (Color online) The net flux as a function of the 
EF frequency. The curves with stars and squares stand 
for results with the EF along x- and z-axes, respectively. 
The dash line denotes the zero-field flux. The window 
amplification within 2-24 THz is exhibited in the inset. 

The extraordinary characteristics of the water trans¬ 
port should uniquely be determined by the behaviors of 
the in-tube water in the applied TEF. We thus count 
the number of hydrogen bonds Nh inside the SWCNT 
in the frequency range as shown in Fig. 3, according to 

O 

the criterion of the oxygen distance less than 3.5 A and 
hydrogen-bond angle < 30°®. It is observed that the Nr 
values are about 3.01, 2.97, 0.81, 0.61, 1.26 and 2.97 at 
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FIG. 3: (Color online) The average occupancy number 
N and hydrogen bond number Nh of water molecules 
inside the SWCNT as a function of the EF frequency in 
the X- or z-direction field, as denoted. 


0.1, 1, 10, 14, 20 and 30 THz, respectively. The min¬ 
imum of the N is about 3.5 at 14 THz. As shown in 
Figs. 2 and 3, the deep valley around 14 THz in the fre¬ 
quency profile of the Nh corresponds to the giant peak 
of the flux and the moderate dip of the number of wa¬ 
ter molecules in the SWCNT, exhibiting a characteristic 
resonant phenomenon. In the whole frequency region in¬ 
terested in this work, we see that the variation of the 
flux takes place coherently with those of the numbers of 
hydrogen bonds and water molecules inside the SWCNT. 
The dramatic drop of the hydrogen-bond number is a re¬ 
sult of the breaking of hydrogen bonds, while the latter 
drives the sudden rise of the net flux. 

Under the equilibrium condition, the in-tube water 
transports in correlated bursts of the hydrogen-bonded 
filament the net flux increases linearly with the pres¬ 
sure gradient. The water transports at equilibrium can 
also be demonstrated by the collective diffusion model®. 
Here, the dramatic rise of the net flux occurs with the 
rupture of hydrogen bonds due to the energy transfer 
from the TEF (see Supporting Material®^). The large 
kinetic energy of water molecules acquired then leads to 
fast shuttling of the water molecules from one end to the 
other and correspondingly the dramatic rise of the flux. 
In this case, the energy transfer from the TEF spoils the 
spontaneous diffusion due to the thermal fluctuation at 
equilibrium, and we cannot follow simply the collective 
diffusion model under the non-equilibrium condition. In 
the frequency range without significant breaking of hy¬ 
drogen bonds, the fluxes in the presence of the pressure 
gradient are consistent with those given by the approach 
by Zhu et al®. In the case of dramatic rupture of hy¬ 
drogen bonds, the relation between the net flux and the 
pressure gradient is not far off linear in a large domain 
of the pressure gradient, see Ref.®^, though in this case 
the equilibrium diffusion does not apply at null pressure. 

Now, we need to reveal the underlying mechanism for 
the fracture of hydrogen bonds. Undoubtedly, the EF 
can’t generate directly the translational motion of charge 
neutral molecules, but can drive the vortical motion of 


FIG. 4: (Color online) The water molecule revolves 
along (a), the axis vertical to the water molecule plane 
through the center of mass, (b), the symmetry axis, (c), 
the axis through the center of mass but parallel to the 
connecting line of two hydrogen atoms, and (d), the 
axis on the hydrogen bond, respectively. 


polar water molecules to acquire the energy. Shown in 
Fig. 4 are four representative vortical modes with differ¬ 
ent rotational axes. The vortical motion can be trans¬ 
mitted to the wobbling of the molecular center of mass 
in the vertical section due to the hydrogen-bond lock¬ 
ing together with the thermodynamical motion. As a 
result, this provides a pathway to vibrational resonances 
to shake off the hydrogen bonding, while the resonant 
frequency range is estimated around 4-14 THz, accord¬ 
ing to the equality of the energies of the water molecule 
binding and the harmonic oscillator^®. Besides the vi¬ 
brational resonances, the TEF can generate direct rota¬ 
tional resonances. Due to the hydrogen bonding of water 
molecules, all the vortical motions should be, in princi¬ 
ple, the reciprocating motions rather than the strict rota¬ 
tion. For convenience, we make somehow discretionary 
identifications for vortical modes in Fig. 4: 4(a), (b), 

and (c) are three independent rotational modes whose 
rotational axes are the principal axes, and 4(d) the li- 
brational mode for the hydrogen bonding. Now, we use 
the rotor model to estimate the resonant frequency for 
rotational or librational motions: Eh = luP' j2 where Eh 
is the binding energy of each water molecule, and I is 
the moment of inertia of the water molecule. If / is the 
total moment of inertia, the energy of a rotor can be re¬ 
expressed by the components on the three principal axes. 
Using the value of Eh = 16 kcal/moP, the resonant fre¬ 
quencies of the water molecule are estimated to be about 
14, 17, 24 and 20 THz, corresponding to the modes in 
Fig. 4(a), (b), (c), and (d), respectively. Besides four 
specific cases, the axes can tilt for the thermal motion. 
Consequently, the change of the moment of inertia results 
in the continuation to neighboring frequency domains of 
the characteristic modes above. 

To evidence these resonant mechanisms, we first do a 
test by extending the TEF to the reservoirs at both sides. 
We find nearly a five-fold enhancement of the net flux in 
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the frequency range of 4-24 THz^ while beyond this fre¬ 
quency range the flux drops off rapidly. The reason for 
this to occur is that the water molecules liberated from 
the hydrogen bonding in the reservoir can break into the 
nanotube rather freely. In addition, the flux enhance¬ 
ment is almost independent of the TEF orientation in 
this case. As the TEF over the reservoirs withdraws, 
some distinct behaviors are exhibited for the single-filed 
water chain in the SWCNT. We see from Fig. 3 that the 
breaking of hydrogen bonds is not efficient below 8 THz 
because of the SWCNT suppression of the needed large- 
amplitude vibration. We also see that the breaking of 
hydrogen bonds through the rotational mode of Fig. 4(c) 
is of low efficiency in that the coherent flip does take 
the priority for the constraint of the SWCNT. In these 
cases, the suppression of the flux occurs correspondingly, 
as shown in Fig. 2. In the a;-direction TEF, an appar¬ 
ent resonant phenomenon, characteristic of a giant peak, 
thus arises eventually. We can see that the resonant fre¬ 
quencies given by the rotor model are consistent with the 
large fluxes and small hydrogen bond numbers shown in 
Figs. 2 and 3. For the ^-direction TEF, we see that the 
rotational mode of Fig. 4 (b) is out of operation. Indeed, 
all the vortical modes in the z-direction TEF are largely 
suppressed by the forced alignment along the z-axis. As 
a result, the giant peak crashes down, as shown in Fig. 2. 

In the low-frequency range roughly below I THz, the 
difference from the zero-field flux can be attributed to 
the alignment of the polar molecules along the EF. The 
change of the alignment direction in the low-frequency 
TEF yields coherently the flips of the water chain, re¬ 
sulting in the moderate enhancement of the flux. With 
the z-direction TEF, the polar water may be aligned 
along the z-axis, the wobbling and random motions in 
the transverse section are largely suppressed with the 
enhancement of the ordered permeation. This interprets 
that the net fluxes remain a little larger than those of the 
x-direction TEF at low frequencies. At higher frequen¬ 
cies (/ > 30 THz), being far away from the resonant 
range, the phenomenon evolves trivially, being close to 
the zero-field results. 

To comprehend the transport behaviors in the TEF 
in more details, we now examine the average dipole ori¬ 
entations of the water chain. Fig. 5 shows the prob¬ 
ability distributions of the average angle (0) of water 
molecules inside the SWCNT for different frequencies, 
where 6 represents the angle between the symmetry axis 
of the water molecule [see Fig. 4(b)] and the z-axis. The 
average is performed over all the water molecules inside 
the SWCNT in last 40 ns of the simulation time. As 
shown in Fig. 5, only for the zero-field case exists the 
single peak with the peak value around 30°, in accord 
with previous studies^®’^^. Note that the symmetric ori¬ 
entation distribution around 30° and 150° at the zero 
field can be anticipated under small pressure gradients, 
because in this case low speeds of water molecules can 
yield sufficient time to produce significant flips in the 
nanotube. With the frequency at 0.1, 1, 24 and 40 THz 
appears remarkably the double-peak structure, serving 
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FIG. 5: (Color online). Probability distributions of the 
dipole orientation of water molecules inside the tube for 
different EF frequencies. 


as the signature of the frequent flip of water molecules. 
The flip may occur for some perturbation of the form 
like in a sufficient time span. The noticeable fact is 
that the double-peak structure becomes flattened for the 
frequencies within the top region of the resonant peak 
in Fig. 2, say, 10 , 14 and 20 THz, as shown in Fig. 5. 
This provides a specific evidence for the hydrogen bond 
fracture in the resonant zone. It is the release of water 
from the hydrogen bonding that smears out the specific 
orientations peaked around 30° and 150°. Note that the 
less flattening for the 20 THz curve is also consistent 
with the farther distance from the peak value as shown 
in Figs. 2 and 3. 

In addition, the resonance can be well visualized in a 
dynamic evolution of hydrogen bonds for water molecules 
confined in the SWCNT. We observe that the average 
numbers of hydrogen bonds at 0.1 and 30 THz remain 
almost unchanged from 0 to 10 ps. However, we observe 
that in the frequency range of 10-20 THz the hydrogen- 
bond number starts to reduce explicitly within 1 ps. 
These exemplifications indicate that resonant response 
of the water chain to the TEF is rather prompt once the 
strong torque is exerted on the water molecule. 

It is worthwhile discussing the dependence of the water 
transport on the TEF amplitude. We find that the aver¬ 
age hydrogen bond number and the net flux do not un¬ 
dergo significant differences for Eq = 1 and a few V/nm. 
However, for the Eq below 1 V/nm, the TEF effect on 
the water transport fades away quickly with the decrease 
of the TEF amplitude. For instance, at Eq = 0.1 V/nm, 
the average hydrogen bond number and the net flux are 
similar to their zero-field results, respectively. The rea¬ 
son of this phenomenon lies in the fact that the small- 
amplitude TEF can not produce sufficient torque to win 
the competition with the thermal fluctuations of the wa¬ 
ter molecules. It is also significant to explore the resonant 
effects on the water transport in much longer SWCNTs 
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that is easier to connect with the experiments. In the 
SWCNT with the length of 10 nm, we have observed the 
similar rupture of hydrogen bonds. The flux enhance¬ 
ment due to the rupture of hydrogen bonds is also about 
3-fold at the same pressure gradient as that for the short 
SWCNT. Moreover, we have studied the water transport 
properties at much smaller pressure gradients. It is ob¬ 
served that the resonant phenomenon is not changed by 
the magnitude of pressure gradient. We perform, as an 
example, the simulation with the frequencies 10 and 14 
THz. Compared to net fluxes at non-resonant frequen¬ 
cies, the significant flux enhancement remains at lower 
pressures. For details, see the Supporting Material^^. 

We notice that the TEF should be regarded as the ef¬ 
fective field of the THz radiation penetrating the tube 
because of the partial screening by the electron clouds 
laying on the SWCNT. It is of interest to consider later 
on whether the radial breathing mode (RBM) of the 
SWCNT can be efficiently generated by the TEF and how 
RBM-induced vibrations mix with the direct TEF-driven 
modes^®. In addition, we would mention that simulations 
with polarizable models for conducting carbon nanotubes 
do not produce appreciable difference in the water trans¬ 
port, compared with those with nonpolarizable models 
for semiconducting carbon nanotubes®^’®®. Finally, in 
order to examine the model dependence, we performed 
some comparative simulations using the SPC/E water 
model®®’®^. Though the magnitude of the flux with the 
TIP3P water model is different from the one obtained 
with the SPC/E model, the tendencies of the net flux, 
the hydrogen-bond and occupancy numbers with both 
models are consistent with each other. Thus, the main 
conclusions that we have drawn do not depend on the 
specific water model. Eor details, see Ref.®^. 


IV. CONCLUSIONS 


In summary, we have demonstrated the effects of 
the vortical resonant mechanisms in the TEF on the 
transport properties and dynamical behaviors of water 
molecules across the SWCNT. The resonances, resulting 
from librational, rotational and rotation-induced trans¬ 
lational modes of polar water molecules, generate signifi¬ 
cant hydrogen-bond fractures and consequently the dra¬ 
matic enhancement of the net flux, with the TEF normal 
to the SWCNT. The resonant phenomenon is well estab¬ 
lished. The change shows up when the TEF turns to 
the symmetry axis of the nanotube: the resonant phe¬ 
nomenon is greatly suppressed by the alignment of water 
molecules along the symmetry axis, which characterizes 
the quasi one-dimensional feature of the SWCNT. Our 
work suggests that the TEF resonant mechanisms are 
instructive to develop high-flux nanoscale devices and ex¬ 
plore possible implications for the biological water chan¬ 
nelling. 
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